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ABSTRACT 

We explore a set of single-point chemical models to study the fundamental chemical aspects of episodic 
accretion in low- mass embedded protostars. Our goal is twofold: (1) to understand how the repeated 
heating and cooling of the envelope affects the abundances of CO and related species; and (2) to 
identify chemical tracers that can be used as a novel probe of the timescales and other physical aspects 
of episodic accretion. We develop a set of single-point models that serve as a general prescription for 
how the chemical composition of a protostellar envelope is altered by episodic accretion. The main 
effect of each accretion burst is to drive CO ice off the grains in part of the envelope. The duration of 
the subsequent quiescent stage (before the next burst hits) is similar to or shorter than the freeze-out 
timescale of CO, allowing the chemical effects of a burst to linger long after the burst has ended. We 
predict that the resulting excess of gas-phase CO can be observed with single-dish or interferometer 
facilities as evidence of an accretion burst in the past 10 3 -10 4 yr. 

Keywords: stars: formation — stars: protostars — circumstellar matter — accretion, accretion disks 
— astrochemistry 



1. INTRODUCTION 

The wide spread in bolometric luminosities ob- 
served among low-mass embedded protostars is at odds 
with a constant mass accretion r ate onto the central 
source (iKenvon et al.1 [1990L fl99i IDunham et al.1 [200l 
lEnoch et all 120091: lEvans et al.ll2009D . Several explana- 



tions have been put forward for this "luminosity prob- 
lem," the most common of which is episodic accre- 
tion: protostars accrete most of their mass in short 
bursts, while spending mo st of their time in a low- 
luminosity quiescent state (jKenvon et al.l 119901 ) . The 
hypothesis of episodic accretion is supported by a 
suite of models (lYoung fe Ev ans 200 51: iVorobvov fc Basul 
20051 120101: iZhu et al.1 120091 120101: IDunham et al.ll2010t 



Dunham fc Vorobvovll2012[ ). but many key questions re- 
main unanswered. How intense are the accretion bursts? 
How many bursts does a given protostar experience, and 
on what timescales do they occur? How does episodic 
accretion affect the chemical evolution from a molecu- 
lar cloud to a circumstellar disk? This Letter focuses 
on the question of episodic accretion chemistry and ex- 
plores how one might use certain chemical tracers to help 
answer the questions outlined above. 

The lack of constraints on most aspects of episodic 
accretion in embedded Class and I protostars comes 
from the difficulty of observing the accretion bursts di- 
rectly. In more evolved Class II sources, ac cretion bursts 
show up as FU Ori and E X Ori events (|Herbigl fl977l : 
lHartmann fc Kenvonl [T996). Luminosity flares of fac- 
tors of 2-1 have been observed for a few late Class 
I sou rces (jKospal et al.l 120071 ; iCaratti o Garatti et al.l 
2011). Other evidence for variable accretion in em- 
bedded s ources is all indire ct, such as the luminosity 
problem (jKenvon et al.lll990l ) o r the presence of periodic 
shocks along protostellar jets (jReipurthl 119891 ) . Addi- 
tional probes are needed, and chemical signatures may 
be one. 



To date, only two studies have addressed the topic 
of episodic accretion chemistry in embedded protostars. 
Lec (2007) showed that repeated burst-quiescent cycles 
lead to repeated freeze-out-evaporation cycles of CO and 
other species. The freeze-out timescale is longer than the 
duration of a burst, so the chemical effects of the burst 
can remain visi ble through part o r all of the subsequent 
quiescent stage. iKim et all (|201lD hypothesized that ev- 
ery time CO freezes out, a fraction is converted into CO2 
ice. They found a better fit to the 15.2 /zm CO2 ice fea- 
ture in the low-luminosity source CB130-1 - presumed to 
be in a quiescent stage - if 80% of CO is converted into 
C02- It is unclear, however, why the conversion would 
have to happen in between accretion bursts rather than 
b efore the o nset of collapse. 

Eei (|2007D and IKim et al.l (|2011l) used a full envelope 
model with fixed timescales for the burst and quiescent 
stages. In this Letter, we explore episodic accretion 
chemistry at an even more fundamental level by way of 
single-point models at various densities and timescales 
(Sects. [5] and [3]). We identify the key chemical pro- 
cesses affected by episodic accretion and illustrate how 
the physical and chemical timescales play against each 
other. We also offer some suggestions on how the chem- 
ical effects of episodic accretion may be observed with 
single-dish and interferometer facilities (Sect. [4j). 

2. THEORY AND MODEL DESCRIPTION 

The fundamental aspects of the effects of episodic ac- 
cretion on the chemistry of embedded protostars are best 
explored with a simple parametric model. Every accre- 
tion burst is accompanied by an increase in the proto- 
stellar luminosity, which raises the dust and gas temper- 
atures throughout the envelope. Once the burst ends, 
the temperatures return to the original quiescent values. 
The chemical response to the temperature changes in the 
envelope depends on the duration of the bursts and qui- 
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escent periods and on the timescales of the dominant 
chemical reactions. The heating and cooli ng timescales 
for the gas and dust are less than a day ([Draind [20031) 
and do not feature into the discussion. 

At typical envelope temperatures of at most 100 K, the 
chemistry following a decrease in temperature is domi- 
nated by the freeze-out of some gas-phase species onto 
the cold dust. Conversely, a heating event results in 
thermal evaporation of some ices. The evaporation rate 
depends exponentially on the dust temperature and is 
independent of the gas density; it is effectively instanta- 
neous (r cv < 1 yr) once the temperature exceeds a given 
species' evaporation temperature. The freeze-out rate 
is gove rned by the collision r ate between molecules and 
grains (|Charnlev et al.l [20011 ). The freeze-out timescale 
of CO and N 2 onto 0.1 /xm grains is 

, /lOK 10 6 cm- 3 

7fc = l x 10 4 y M / — — — , 1 

V T g n(H 2 ) 

for a sticking coefficient of unity (Bissch op et al.l l2006) . 

The chemical timescales have to be balanced against 
the duration of the accretion bursts (rb) and of the quies- 
cent p eriods betw e en sub seque nt bursts (r a ) . In t h e mod - 
cls of iZhu et al.l (pOlOl ) and IVorobvov fc Basul (pOlOh . 
bursts typically occur every 10 3 -10 4 yr. The duration 
of each burst is at most a tenth of the duration of the 
quiescent phase. For some initial conditions, the models 
predict no substantial bursts at all. 

Having identified the relevant timescales, the next 
question is how much the envelope is heated during 
a burst. At early times, appropriate for embedded 
protostars, the lu minosity is do minated by accr etion 
ilMvers et al.|[l998h The models of lZhu et all tpQIOh and 
IVorobvov fe Basul (|2010l ) have mass accretion rates of 
M q = 10~ 9 -10 -7 Mq yr -1 in the quiescient phases and 

of Mb = 10 -6 -10 -4 Mq yr -1 during accretion bursts. 
For a stellar mass of 0.3 Mq and a stellar radius of 2 
i?0, the corresponding accretion luminosities (L a cc ~ 
GM*M/R*) are 0.005-0.5 and 5-500 L & , respectively. 
We use the continuum radiat ive transfer code RADMC 
(jDullemond fc Dominik 2004|) to compute dust temper- 
ature profiles for an envelope exposed to central lumi- 
nosities between 0.005 and 500 Lq. Our model envelope 
extends from 20 to 10 4 AU and has a mass of 1 Mq, dis- 
tributed along a power-law density profile with an expo- 
nent of— 1.5. The gas-to-dust mas s ratio is set to 100 and 
the du st opacities are taken from [Osscnkopf fc Henninel 
(1994) for thin mantles at 10 6 cm" 3 (OH5). The goal 
is to get a rough idea of how much the temperature in- 
creases during a burst, so the exact choice of envelope 
parameters and dust properties is not important. 

Figured] shows the temperatures at the radii where the 
density is 10 5 , 10 6 , and 10 7 cm -3 . A lower limit of 10 
K is applied in the low-luminosity cases to account for 
external heat ing by cosmic ra ys and the interstellar ra- 
diation field ([Galli et al.ll2002f) . An accretion burst with 

an intensity of Jb = Mb/Af q = 100 or 1000 increases the 
temperature by a factor of 2.5 or 4, respectively. This 
relationship holds regardless of the density power-law ex- 
ponent and other envelope parameters. 

We explore the chemical aspects of episodic accretion 
in a set of single-point models at densities of 10 5 , 10 6 , 
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Figure 1. Dust temperatures in a spherical envelope as function 
of the stellar luminosity and the mass accretion rate. The three 
curves correspond to the radii where the density reaches 10 s , 10 6 , 
and 10 7 cm _i . 

and 10 7 cm -3 . For each de nsity, the dust an d gas tem- 
peratures (assumed coupled; [Galli et aT]|200 2) are cycled 
between 10 and 30 K, corresponding to an accretion burst 
with an intensity of I\, = 240. The temperatures of 10 
and 30 K fall below and above th e CO and N 2 eyapora - 
tion temperatures of about 20 K ([Bisschop et al.l [20061. 
so we expect the overall chemical effects to be dominated 
by their freeze-out-evaporation cycles. 

We adopt the UDfA06 chemical network 
(|Woodall et al.l 120071) . expanded with freeze-out and 
thermal evaporation of all neutral species. The binding 
energies of CO, N 2 , H 2 0, NH 3 , a nd CQ 2 are set to 
855, 800, 5773, 2790, and 2300 K (iFraser et al.1 [20(51 
Bolina fc Brownll2005l: [Bisschop et al.ll2006l: iNoble et al.l 
20121 ). We follow the chemistry for five burst-quiescent 
cycles with Tb = 100 yr at 30 K (burst) and r q = 10 3 
or 10 4 yr at 10 K (quiescent). The computation 
begins with typical dense pre-stellar core conditions as 
constrained by observations and models ([Caselli et al.l 
[19991: IGibb et al.|[200l IMaret et al.|[2006[ ): hydrogen in 
H 2 ; carbon in CO (80%), CO ice (10%), and C0 2 ice 
(10%); remaining oxygen in H 2 ice; and nitrogen in 
atomic N (60%), N 2 (20%), N 2 ice (10%), and NH 3 ice 
(10%). Also included are He, S, Si, Fe, N a, Mg, P, and 
CI. E lemental abundances are taken from Aikawa et al.l 
( 2008) and the cosmic-ray ionization rate of H 2 is set to 
5 x 10" 17 s" 1 (|Dalgarnoll2006l) . 

3. RESULTS 

The temporal variation of the gas-phase abundances 
of CO, N 2 , HCO+, N 2 H+, and H+ is plotted in Fig. 
[H The six panels show how the chemistry is controlled 
by the interplay of the dominant physical and chemi- 
cal timescales. In panel f, the freeze-out timescale of 
CO and N 2 is the same as the duration of the quiescent 
phase (rf r = r q = 10 4 yr). Both species deplete by sev- 
eral orders of magnitude before the next burst hits. If 
we decrease the density by a factor of 10 or reduce the 
spacing between bursts by the same factor (panels c and 
e), freeze-out affects only half of the available CO and 
N 2 . At even lower densities (panels a, b, and d), the 
freeze-out timescale is 100 or 1000 times longer than the 
duration of the quiescent phase, so CO and N 2 remain 
fully in the gas phase. 
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Figure 2. Abundances of CO, N2, HC0 + , N2H+, and Hg as function of time through five episodic accretion cycles in single-point models 
with bursts every 10 3 or 10 4 yr (top to bottom) and with densities of 10 5 , 10 6 , or 10 7 cm -3 (left to right). The dotted lines represent H^~; 
the other lines are as indicated in panels a and d. 



The chemistry of HCO+, N 2 H+, and HJ is closely tied 
to that of CO and N2. H^" is destroyed mainly by reac- 
tions with CO and N2, which act as the primary for- 
mation channels for HCO + and N 2 H + . The primary 
destruction channels differ for HCO + and N 2 H + : for 
the former it is dissociative recombination with elec- 
trons, producing CO and H2, and for the latter it is 
proton exchange with CO, producing HCO + and N2. 
All these reactions are fast, with timescales of about 
0.1 yr [10 6 cm _3 /n(H 2 )] at the beginning of each quies- 
cent phase. The episodic accretion chemistry of HCO + , 
N2H + , and is therefore controlled by the freeze-out 
timescales of CO and N2. 

HCO + , N2H + , and have the most complex abun- 
dance profiles in panel f of Fig. [5J where the dominant 
chemical timescale (freeze-out) is equal to the spacing 
between subsequent bursts. In each cold stage, H^" be- 
comes more abundant as CO and N2 freeze out. The 
reverse happens during each burst: the rapid evapora- 
tion of CO and N2 results in a rapid loss of Hjj\ The 
abundance of HCO + remains constant between bursts 
as long as the decrease in CO is balanced by an increase 
in H^". When the abundance levels off, HCO + can 
no longer be formed as efficiently, and its abundance de- 
creases. HCO + is reformed during the bursts from the 
evaporation of CO. The third species, N2H + , shows a 
humpback pattern in every cold phase. This pattern re- 
sults from the depletion of CO early on in the quiescent 
phase (slowing destruction of N2H + ) and the levelling 
off of combined with ongoing depletion of N2 at later 
times (slowing formation of N2H + ). 

If the model parameters are changed so that Tf r = 10r q 
(panels c and e) , only half of the gas-phase CO and N2 



still freezes out in between bursts. The abundance pat- 
terns of H^" and N2H + are now inversely proportional to 
those of CO and N2, and HCO + remains constant except 
for small variations due to the inverse t emperature de - 
pendence of dissociative recombination (Mitchell 1990). 
In panels a, b, and d, the dominant chemical timescales 
are so much longer than the physical timescales that all 
abundances remain constant in time except for the same 
small variations in HCO + . 

The initial conditions have little effect on the above 
results. All CO and N2 ice evaporates in the first burst at 
t = 0, so it does not matter what fraction of each species 
starts as gas or as ice. The ratios of CO versus CO2 and 
of N 2 versus NH3 affect the absolute abundances of CO 
and N2 , but leave the nature of the episodic abundance 
variations in Fig. [^unchanged. 

Although these single-point models are only a crude 
analog to an actual protostellar envelope, they do high- 
light two fundamental aspects of episodic accretion chem- 
istry: timescales matter, and the "chemical temperature" 
tends to be higher than the dust temperature. The first 
point is obvious from Fig. [5] the abundances only show 
strong variability if the duration between two accretion 
bursts is comparable to or exceeds the dominant chemical 
timescale, which in this case is the freeze-out timescale 
of CO and N2. The second point is addressed in more 
detail in the next section. 

4. DISCUSSION 

4.1. The chemical temperature and implications for 
observations 

Whenever the stellar luminosity changes, the temper- 
ature of the gas and dust throughout th e envelope re - 
sponds on timescales of at most a day (Draine 2003). 
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Figure 3. Abundances of CO, CO2, CO2 ice, and N2H+ as function of time through five episodic accretion cycles in single-point models 
with bursts every 10 4 yr and a density of 10 7 cm -3 . The dotted lines represent CO2 ice; the other lines are as indicated in panel a. The 
panels differ in the fraction of CO ice converted into CO2 ice in each quiescent stage. 



The chemistry responds equally fast to heating events, 
but it takes at least 10 3 yr to adjust to a drop in tem- 
perature. This can lead to a prolonged situation where 
the dust and gas temperatures lie below the evaporation 
temperature of a given species, but the species remains 
in the gas phase because it has not had enough time 
yet to freeze out. Obscrvationally, we therefore predict 
a mismatch between the temperature inferred from dust 
continuum emission and the "chemical temperature" in- 
ferred from the molecular abundances. In other words, 
the molecule does not exhibit the depleted ab undance 
profile expected based on the dust temperature. lAikawal 
(2007) employed similar arguments to explain how verti- 
cal mixing in a circumstellar disk can sustain a reservoir 
of gas-phase CO below 20 K. 

As an example, in panels a-c of Fig. [2j each model 
spends 90% of the time at a temperature of 10 K and 
10% at 30 K. Observations of the dust emission would 
recover those fractions correctly. However, CO spends 
most of the time in the gas phase, so observations of the 
CO emission would imply a temperature of 30 K nearly 
100% of the time. In the absence of accretion bursts, 
this dense gas should be dominated by freeze-out. Hence, 
the mismatch in temperatures inferred from the dust and 
from the gas-phase abundances can be used as a novel 
probe for episodic accretion in embedded protostars. 

The situation is more complicated in reality because a 
protostellar envelope spans a range of densities and tem- 
peratures: typically ab out 10 4 -10 9 cm -3 and 10-250 K 
(jJeirgensen et alJl2002f >. In a given accretion burst cycle, 
only a limited part of the envelope is heated from below 
to above the CO evaporation temperature. Nevertheless, 
our fundamental conclusion about relative timescales still 
applies: if the duration between bursts is shorter than the 
freeze-out timescale in the affec ted part o f the envelope, 
CO remains in the gas phase (jLed 12007). Hence, if an 
envelope is observed between bursts, we expect to find 
an excess of gas-phase CO relative to the inferred dust 
temperature profile. 

How could one probe such an excess? The CO abun- 
dance profile is commonly constrained by single-dish ob- 
servations of pure rotational emission lines from th e op- 
tically thin C ls O isotopolog (pOrgensen et al.ll2002t ). By 
combining spectra from a series of rotational lines with 
dust continuum data, it is possible in principle to infer 
the temperature at which CO freezes out. The excess CO 



gas exists at temperatures below the actual freeze-out 
temperature, so the temperature inferred from the ob- 
servations should be lower t han the value of abou t 20 K 
measured in the laboratory ( Bisschop et al.ll2006f ). How- 
ever, even with access to seven of the ten C 18 lines from 
J = 1-0 to 10-9, the freeze-out temp erature derived from 
observations is poorly constrained (jYildiz et al.l 120101) . 
Single-dish observations of CO isotopologs in any one 
particular source are therefore unlikely to be a viable 
probe of episodic accretion chemistry. 

We suggest two alternative observational probes of ex- 
cess CO due to episodic accretion. The first is to take 
a large sample of prot ostars, such as th at from the c2d 
Spitzer legacy project (jEvans et al.ll2009l ) , and search for 
correlations between low- J C 18 line intensities (directly 
tracing the column of cold CO) and systemic properties 
such as bolometric temperature, bolometric luminosity, 
and envelope mass. The second method is to directly 
measure the size of the C 18 emitting region in a sample 
of protostars using interferometers such as the SMA or 
ALMA. If there was a recent burst, the size of the CO 
emitting region should extend into dust with tempera- 
tures below the freeze-out temperature. Either method 
would need to be coupled to a parametric study of line 
radiative transfer models to understand in more detail 
how the low-J C 18 emission - spatially resolved and 
unresolved - depends on variations in Tboi, ^bol, and 
M cnv , as well as on variations in the episodic accretion 
timescales. 

4.2. Grain-surface chemistry of CO 

iKim et all (|2011[ ) suggested that a fraction of the 
CO ice is converted into CO2 ice in every quiescent 
stage. CO2 has an evaporation temperature of 35-40 K 
( Noble et al. 2 0l2|) . If the conversion takes place in part 
of the envelope where the quiescent temperature lies be- 
low the CO evaporation temperature, it is unlikely to 
get warm enough during a burst for CO2 to evaporate 
(Fig. Q]). The presence of a pure CO2 ice feature at 15.2 
/im in low-luminosity sourc es may th erefore indicate an 
accretion burst in the past (|Leell2007l) . 

We reran the model for n(H2) = 10 7 cm~ 3 and r q = 
10 4 yr (Fig. |2J), converting either 20, 50, or 80% of CO 
ice into CO2 ice in each quiescent stage. The resulting 
abundances of CO, CO2, CO2 ice, and N2H + are plot- 
ted in Fig. [3l Regardless of the conversion fraction, the 
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abundance of CO2 ice rapidly goes up to about 10~ 4 . 
The gas-phase abundance fluctuates between 3 x 10" 11 
and 1 x 10" 10 in response to the abundance variations of 
H3" (Fig. [2]). N2H + gradually becomes more abundant in 
each burst-quiescent cycle as its main destroyer, CO, is 
lost. 

If a substantial amount of CO is indeed converted into 
CO2, we would no long er expect to see an excess of 
gas-phase CO (Sect. 14. ip . but instead we would expect 
a higher abundance of N2H+ between bursts. Conver- 
sion of CO ice into other species, such as H2CO and 
CH3OH, would have the same effect. At lower densities 
or a shorter quiescent stage, the amount of CO ice built 
up in each quiescent phase is smaller than for the case 
explored in Fig. [3l and the effects of grain-surface con- 
version of CO into CO2 would be subdued. A grid of full 
envelope models is needed to explore quantitatively how 
much pure CO2 ice is formed as function of various pa- 
rameters. As it stands, our single-point models suggest 
that the abundances of CO2 ice and gas-phase CO and 
N 2 H+ together are a novel probe to elucidate episodic 
accretion in low-mass embedded protostars. 

5. CONCLUSIONS 

Low-mass protostars appear to accrete most of their 
mass in a series of bursts interspersed with relatively 
long phases of very little accretion. If this hypothesis 
of episodic accretion is indeed true, the protostellar en- 
velope is exposed to multiple heating and cooling events 
over the course of its existence. This Letter presents 
a fundamental theoretical investigation into how this 
would affect the chemical composition and evolution of 
the gas and dust in the envelope. 

Based on simple single-point models, we conclude 
that the effects of episodic accretion chemistry depend 
strongly on the relative physical and chemical timescales. 
The abundances of common species like CO, N2, HCO + , 
and N2H + show strong time variability if the duration 
between subsequent accretion bursts is comparable to 
or longer than the dominant chemical timescale. For 
typical cold envelope conditions, the dominant chemical 
timescale is the freeze-out timescale of CO and N2. Re- 
gardless of the exact timescales, the chemical effects of an 
accretion burst should remain visible for 10 3 -10 4 yr af- 
ter the burst has passed. In particular, we predict excess 
amounts of CO gas to be observable with either single- 
dish or interferometer facilities as chemical signposts of 
episodic accretion. 

This work was supported by the National Science 
Foundation under grant 1008800. 
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